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Abstract. For smooth domains, Liu et al. (Comm. Pure Appl. Math. 60: 1443- 
1487, 2007) used optimal estimates for the commutator of the Laplacian and the 
Leray projection operator to establish well-posedness of an extended Navier-Stokes 
■ dynamics. In their work, the pressure is not determined by incompressibility, but 

^v^j \ rather by a certain formula involving the Laplace-Leray commutator. A key estimate 

^ ■ of Liu et al. controls the commutator strictly by the Laplacian in L 2 norm at leading 

| order. In this paper we show that this strict control fails in a large family of bounded 

planar domains with corners. However, when the domain is an infinite cone, we find 
\ that strict control may be recovered in certain power-law weighted norms. 

< 

In this paper, we study estimates for [A, P] = AP — PA, the commutator of the 
Laplacian and the Leray projection operator, in planar domains with corners. In 
a bounded domain Q C M , the Leray projection operator P is defined as follows: 
Given any a G L 2 (Q,R N ), there exists a unique q G H 1 ^) with f n q = and such 
that Pa := a + Vq satisfies 

(1.1) = (Pa, V0) = (a + Vq, V(j)) 



in 

m 



1. Introduction 



for all G H 1 ^). In [3], Liu et al. proved the following L 2 -estimate for the commu- 
tator of the Leray projection operator and the Laplacian. 

Theorem 1.1. Let Q be a connected, bounded domain in M. N , N > 2, with C 3 
boundary. For any f3 > |, there exists C > such that for all vector fields u G 

>< . //*n// ( i(a E 

. <£; (1.2) / \[A,P]u\ 2 </3 I |Au| 2 + C7 / |Vu| 2 



Jn Jn Jn 

Theorem II. II has significant applications to the Navier-Stokes equations. We recall 
that on a bounded domain Q in M> N for N > 2, the Navier-Stokes equations modeling 
incompressible viscous fluid flow with no-slip boundary conditions are given by 

{d t u + u ■ Vu + Vp = vAu + / 
V-u = 
u\ r = 0, 

where T = dQ, u denotes the velocity of the fluid, p denotes the pressure, and v 
represents the viscosity. In the authors consider strong solutions to (NS) with 
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constant v > 0, and show that the pressure satisfies 

(1.3) Vp= (I - P)(f -u-Vu) + u[A,P]u. 

For such solutions they prove the unconditional stability and convergence of a simple 
time discretization scheme which decouples the updates of velocity and pressure. The 
decoupling of these variables is significant in that it eliminates the need for an inf-sup 
condition which is often necessary to prove the stability in finite-element schemes. 
A critical ingredient in the proof of stability in [3] is that by invoking Theorem 11.11 
with < 1, one can strictly control the pressure gradient by the viscosity term plus 
lower-order terms. As a result, Liu et al. establish the well-posedness of an extended 
Navier-Stokes dynamics in which the pressure p is always determined by the formula 
( 11. 3j) and the zero-divergence condition is dropped in general. We refer the reader to 
[3] for further details and discussion. 

Theorem 11.11 assumes that the boundary r of Q is C 3 . One would like to weaken 
this assumption to allow, for example, sharp corners on T. In this paper, we show 
that such an improvement is not possible. We let /C CT denote an infinite cone centered 
at the origin, taking the form 

(1.4) K„ = {(xi,x 2 ) G M 2 : < r < oo, < 9 < a}, 

where r and 6 denote the polar coordinates of (xi,x 2 ) and a G (0, 2tt). We consider 
bounded domains 0c8 2 satisfying the following property: there is a neighborhood 
U of such that U n Cl — U PI K, a for some rotated translate Cl = R(Q — xq) of Q and 
for some a ^ n. In this case we call Q a bounded domain with a straight corner. We 
claim that Theorem 11.11 fails on any such domain. 

Theorem 1.2. Let Q in IR 2 be a bounded domain with a straight corner. Then for 
every (3 < 1 and for every C G IR, there is a vector field u G H 2 r\Hl(Vt,M?) satisfying 

(1.5) / \[A,P]u\ 2 >/3 [ \Au\ 2 + C [ \Vu\ 2 . 

Jn Jq Jq 

One may suspect that the reason < 1 is not possible in general has something to 
do with the lack of H 2 regularity for the Stokes operator in domains with reentrant 
corners. One known way of dealing with this situation involves using weighted Sobolev 
spaces. In a recent paper of Rostamian and Soane [I], the authors reformulate the 
time discretization scheme of [3] in non-convex polygonal domains using such weighted 
spaces. While the authors do not prove convergence of their scheme, they do give 
numerical evidence suggesting that this scheme converges to the correct solution. 

We are motivated by jl] and elliptic regularity theory with weights [2] to allow for 
corners on T and look to prove an optimal estimate similar to (II. 2p in a weighted 
L 2 -space. For the most part, we study conical domains of the form in (ll.4p . The 
weighted spaces considered in [2] are defined as follows. 

Definition 1.3. For an integer I > and a real number a, we define the space 
V^ Q (/C CT ) to be the closure of C , ^°(/C cr \{0}) with respect to the (scale-invariant) norm 

(1.6) ll/lk, a = (f K ^r 2 ^' l+ ^\D^J\ 2 rdrde^ < oo. 

We refer the reader to [2j for a more thorough discussion of weighted Sobolev spaces 
in an infinite cone. 
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Before we state the main theorem, we must define the Leray projection operator 
on unbounded domains. This definition differs from that given in (II. ip . because if O 
is unbounded, then Vi? 1 (fi) is not closed in L 2 (Q). To remedy this, we fix a bounded 
domain BcflC R , and we define the space 



(1.7) Y={qe Ll c (Q) : Vg G L 2 (fi, R N ) and / q = 



B 



Then Y is a Hilbert space with norm ||g||y- = f n |Vg| 2 , and the space VY is closed 
in L 2 (f2,M iV ). We define the Leray projection operator P as in (II. ip . except that we 
assume g is in Y(Q) instead of if 1 (f2). Further discussion of the Leray projection 
operator on unbounded domains can be found in [5]. 

We remark that if Q is Lipschitz, C£°(f2) is dense in Y. The proof of this fact is 
similar to the proof for Q = indicated in |3j, based on the case Q = M. N treated 
in [51 Lemma 2.5.4]. 

We are now prepared to state the main theorem. 



Theorem 1.4. Suppose a G (0, 2tt) and let )C a be an infinite planar cone as in ( T^J 
Let q^I. Then the following estimate holds for all u G C^lJCa^O}, M 2 ): 



1.8) I r 2a \[A,P}u\ 2 rdrd9 < (3 a , a I r 2a \Au\ 2 rdr 

/3 a>a = sup max I j3 +jk , /L fe [ , 



d9, 



where 



with 

(i. 9 ) = sr a \d - e — )R ' (i - a + ± 



2k 2 (l-ay ' \(l± e -( k - ia y)(l - e- 2 ( k+ia -^) 

Moreover, f5 Cj0t is the smallest constant satisfying < \1.8t\ for every u G C^°(/C (J \{0}, M 2 ). 

We will prove Theorem ll.4l in Sections [2] and |3j In Section [5j we show that Theorem 
11.41 implies Theorem 11.21 

The expressions in (11.91) are sufficiently complicated that it is difficult to charac- 
terize exactly when f3 a ^ a < 1 holds. We will make a few observations, however, and 
provide numerical evidence which suggests that for all a G (0, 2n) except for one value 
o~ = o c ~ 1.43037T, we have {3 a>ol < 1 for a in some interval just to the left or right of 
a = 0. 

First, note that as k — > oo we have P±,k \- For a = we compute that 

- 1 cosh 2 ka — cos 2 a =F cosh ka sin 2 a =F k sin a cos a sinh ka 

(1.10) fi ±k = , 

2 cosh ka — cos^ a 

from which we see that if a = 7r, then j3± t k = |, hence /^o = |. This half-space 
estimate (II. 8p with constant weight was already proved in [3], and explains why the 
condition /3 > \ is essentially optimal in Theorem 11.11 Note that due to the dilation 
invariance of the domain, no lower-order term such as that in (11. 2p should appear in 
the half-space case, since it would scale differently under dilation. 

Whenever ix ^ a G (0,27r) however, we have (3- t Q = 1, /3 +i o — 0. Thus, whenever 
the weight is constant (a = 0) and the cone has a corner (a ^ ir) we conclude that 
the optimal constant /3 CTi o > 1- Our proof of Theorem 11.21 relies on this fact. 
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It is easy to approximate fi a fX numerically. For a number of values of the cone 
angle a, in Figure 1 we plot log 10 a>a vs. a for a G [—1,1]. Spikes appear in many of 
these graphs, providing evidence of singularities where presumably j5 Cj0t = +oo. After 
closer examination, these graphs suggest that: 

• Ax,a = 1 whenever a = and o ^ tt. 

• Pa,a < 1 for small a > when < a < tt or a c < a < 2tt. 

• fla,a < 1 f° r small a < when 7r < cr < er c . 

The number a c ~ 1.43037T satisfying cr c cotcr c —1 = appears to be a critical value 
of o where the minimum of ft a>a occurs at a = 0, and the minimum value is 1. To 
see this, we observe that numerical evidence indicates that for a near the critical 
value and for a near 0, (3 a}0l is achieved at k — 0. We therefore take the limit as k 
approaches of f3 + ^ and /L^, which yields the formulas for /3 +0 and /§_o given in 
( 14. 2 j) and (|4.3|) . Numerical evidence again shows that for a in a neighborhood of the 
critical value and for a near 0, (3 +t o > /3_q. Using Maple to differentiate (3 +t o with 
respect to a, and evaluating the derivative at a = 0, we find that 

(1.11) a /3+,o|a=O = 0-COt<7 - 1. 

These numerical results also suggest that < 1 for convex cones, uniformly in a 
for positive a in a fixed interval. So we may conjecture that for a bounded polygonal 
domain Q that is convex, say, an estimate of the form 

(1.12) / r 2Q |[A,PH 2 < f r 2a ((3\Au\ 2 + C\Vu\ 2 ) 
Jn Jn 

will hold for some /3 < 1 and C independent of u in a suitable space of functions 
vanishing on dfl, provided a is small and positive. Here r = r(x) would be the 
distance from x G Vt to the nearest corner on V. The lower order term on the right 
hand side of (11.121) comes from the definition of the V 2 a norm on a bounded domain 
(see [2]), given by 




We do not include the term IIt""-^! 2 ^ on the right hand side of ( 11.121) . because it can be 
controlled by first order partial derivatives using a Hardy inequality (see [2], Chapter 
7 for details). 

However, we have no proof of f II . 1 2 j) at this time. 

2. Preliminary transform in radius 

From the pressure formula (ll.3p we see that the commutator [A, P]u represents the 
contribution of the viscosity term to the Navier-Stokes pressure gradient. Specifically, 
[A, P]u represents the pressure gradient for the linear Stokes equations with no-slip 
boundary and without forcing. For this reason, as in [3], we refer to the corresponding 
pressure as the Stokes pressure, denoted ps = Ps(u). From (II. ip . when a = u G H 2 (Q) 
with Q unbounded, we have VAg = AVg = VV ■ a and it follows easily (as in [3]) 
that 

[A, P}u = (I - P)(Au - VV • u) = V Ps . 
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FIGURE 1. log 10 (/3 (7)Q ,) vs. a for various o. From left to right, top to 
bottom, cx/tt = 0.2, 0.4, 0.5, 0.65, 0.85, 0.95, 1, 1.05, 1.2, 1.4, 1.6, 1.8. 

We recall from [21 Sec. 2.1] that the Stokes pressure p$ is determined (up to constant) 
as the solution to the boundary value problem 

(2.1) Ap s = in tt, n- Vp s = n • (A - VV-)w on V. 

(The boundary condition holds in H^^iY) due to a standard trace theorem, since 
the vector fields Am — VV • u and Vps are in L 2 (Q, R ) with zero divergence.) 
Letting 



\Vp 



s \\V° 

" V 2,a 



„2» 



|Vp s |V dr d6, 



lAd 



V 2,a 



" 2a \Au\ 2 rdrd6, 



we see that in order to prove Theorem ll.4[ we must determine the smallest constant 
/3 a>a satisfying the inequality I p < (3 ata I u , subject to ( 12.1ft . In this section, we perform 
the first steps in our attempt to find fi a ,u- These steps amount to taking a Mellin 
transform of the problem. We first rewrite I p , I u and (12. ip in terms of the polar 
coordinates {r,8), then change variables using r = e s , which transforms lC a to an 
infinite strip S. Taking a Fourier transform will reduce the problem to a family of 
maximization problems parametrized by a Fourier variable fcGl. 
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2.1. We begin by letting 

— 1 \ jq _ I cos 6 — sin 9 



J \ 1 ) ' 6 \smf) (csfV 
A straightforward calculation shows that 

allowing us to rewrite 7 p as 

I p = [ (\rd r p\ 2 + \dgp\ 2 )r 2a - 1 drd9. 

We change variables by letting r = e s , resulting in a transformation of the domain 
K a to an infinite strip S = {(s,9) G 1R 2 : — oo <s<oo, 0<#< a}. We then let 
q = e as p s and express I p in terms of q. We conclude that 

p poo 

(2.2) I p = (\d s q-aq\ 2 + \d e q\ 2 )dsd6= / I p>k dk, 

J S J -oo 

where k is the Fourier variable corresponding to s, and 

(2.3) I Ptk = [ (\(k + ta)q\ 2 + \d e q\ 2 ) dB. 

Jo 

2.2. To rewrite I u , we first calculate 

Am = V • Vu = (rd r + 2)( r - 1 d r u) + d e {r- 2 d e u) = r~ 2 ((r<9 r ) 2 + d 2 e )u. 
If we let u = re J0 v, we can show that 

(2.4) Au = r- l e je ((rd r + ifv + (d e + J) 2 v) . 

We again change variables to express I u as an integral over S. We let w = e sa v, and 
we find that 

I u = J e 2sa \(d s + l) 2 v + (d e + J) 2 v\ 2 dsdd 

(2-5) „ nQO 

= / \(d s + 1 - a) 2 w + (d e + J) 2 w\ 2 dsd6= / J Ui 

J S J-oo 



oo 



where 

(2.6) I u ,k= [ \(ik + 1 - a) 2 w + (d e + Jfw\ 2 d9. 

2.3. As with I p and I u , we wish to rewrite (12. ip in terms of fc, ^, q, and iy. We perform 
a change of variables and rewrite the first condition of (12. ip as e~ 2s (d 2 + d 2 )p = 0. 
Recalling that q = e sa p, we find that Aq — e~ 2s (2ad s — a 2 )q = 0, so (dj + d 2 s — 2ad s + 
a 2 )q = 0. To rewrite the boundary condition of (12.11) . we observe that the left hand 
side can be rewritten using the equalities (eg, Vp s ) = (e 2 , r _1 (r9 r p s , dep s )) = r~ 1 dep s . 
For the right hand side, we use (12 ,4p . combined with the equality 

_ -i je ( r <9rV • u\ _ _! jg ( rd r ((rd r + 2)v x + d B v 2 )\ 
VV ' M_r 6 \deV-u)- r 6 \de((rd r + 2)v 1 + dev 2 )) 

and the property v = when 8 = and 9 = a, to conclude that 

(eg, Au — VV ■ u) = r _1 ((r(9 r ) 2 + 2rd r )v 2 — ded r v\ 

= -r~ 2 v 2 - d e d r vi = —dgd r V\ 
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for 9 = and 9 = a. We can therefore rewrite the boundary condition in ( 12.1 1) as 
(2.7) d e p s = -dgdffVx. 

Using the equality w = e sa v, we see after a calculation that we can recast ( 12.11) as 
the following boundary value problem on S: 

(d 2 e +d 2 s -2ad s + a 2 )q = in S, 
1 = —dg(d s — a)w\ when 9 = 0, a. 



Finally, taking the Fourier transform of ( 12.81) in s, we have that for each k G R, q 
must solve the boundary value problem 

^ ^ d 2 e q = (k + ia) 2 q for < 9 < a, 

deq = —dgiik — a)vb\ when 9 = 0, a. 

3. Optimization in angle 

In this section, we determine ft aj(X = sup j- subject to (12. 9ft and the no-slip boundary 
condition. First, for fc ^ we suppress the a and a variables and define 

(3.1) p k = sup { : (EH]) holds, and w = for 9 = 0, a 

I Iu,k 



Note that since w is real, we have w(—k,9) = w(k,9) } hence I u ,-k = Iu,k from (12. 6p . 
and similarly I p -k = I P ,k from ( 12.31) . We conclude that (3k is even in k. We define 
(3 Ut0l = sup fc>0 (3k, and we observe that 



/oo /*00 
I P ,k dk < (3 k I u ,k dk < (3 a , a I u . 
-oo J — oo 

We will prove Theorem 3 by computing that (3k = max{/3 +j fe, $- t k} as given by ( 11.91) . 
and by showing that /3 CT>a < /3 CT)Q ,. Since evidently (3 a , a > f3 a ,a, the result will follow. 

3.1. We first rewrite the quantity from ( 12.61) to diagonalize the matrix involved. 
We define 

Then letting — = Vy with y = (2/1, 2/2), and using = V(iA), we rewrite in 
the following way: 

Iu,k = [ \{(tk + l-a) 2 + (d 2 + 2Jd e -l))(Vy)\ 2 d9 
Jo 

= 2/| ((life + 1 - a) 2 + (d 2 + 2iAd e - 1)) y\ 2 d9 
Jo 

(3.3) = 2 f (\L m \ 2 + \L 2 y 2 \ 2 )d9, 

where 







L x = (iib + 1 - a) 2 + d 2 - 2id e - 1, 
L 2 = (ik + l-a) 2 + d 2 + 2id e -l. 
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3.2. We next express the quantity from (I2.3P in terms of the boundary data from 
(j2~lD . From ([ZS]) it is clear that explicitly 

q(k,9) = a + e {k+ia){e - a) + a_e- {k+ia)e , 

[ 1 ' d e q(k, 6) = (k + za)a+e (fc+ia)(9 - ff) - (k + ia)a.e- {k+ia)6 , 

for some complex constants a+ and a_. If we define 

(3.5) u = e - (k+ia)a 
for convenience, we see from ( 13.41) that 

(3.6) $ (M) = /«+ + a_ W) * = 

v y v ' la + u + a_, = 0, 

and 

(3.7) d e q{K9) = \f h l ia \^- a -"\ 6 = a n 

Combining (13. 7p with the equality — zi&x = (Vy)i, we can rewrite the boundary con- 
ditions in (j2.9j) as 

/ 38 n "+ - a-w = + y 2 ), = <r, 

a + u - a_ = d e (y 1 + y 2 ), 9 = 0. 

These equations will be used later to determine a + and a_ from y (note u 2 ^ 1). To 
rewrite I p ^, we apply (12. 9ft and integrate by parts. This gives 



|(<9 S — a)q\ 2 ds d6 = I (k + ia)q(k — ia)q dk d6 
Js 

(dgq(cr)q(a) — dgq(0)q(0)) dk — / dgqdgqdkdO 



s 

oo 



2{ika — a 2 )qq dk d6, 
which, in light of ( 12. 2j) . allows us to write 

(3.9) I p>k = d e q(a)q~(a) - d e q(0)4(0) - f 2(ika - a 2 )qq~ d6. 



In order to write f° 2{ika — a 2 )qq d6 in terms of a + and a_, we use (13. 4p to evaluate 
the dot product and integrate. We conclude that 

£ 2{ika - a 2 )qqdQ = {ika - a 2 ) + (l - e' 2ka ) 

+ (ik - a) {ia_a + + ia + a_) e ^ k+ia ^ (l - e 2iaa ) 



+ (A; + ia)(a__a + + a + a_)(c<; — u). 

Similarly, to compute doq(o)q(cr) — dgq(0)q(0), we use the formulas for q and dgq on 
the boundary given in ( 13.61) and (13. 7p to write 

n d e q(a)^(a) - 3 e g(0)$(0) = ( k + ^){{\a + \ 2 + |a_| 2 )(l - \u\ 2 ) 

+ (a+a:_ + a^a + )(u — oj)}. 
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Plugging ( 13.1 Op and (13. lip into (I3.9p . we discover that 



(*^) (KI 2 + KI 2 ) (i-H ! ). 



3.3. By the results of the previous subsection, (3k is the supremum of the ratio 
I P ,k/I u ,k subject to (13 .8p and the no-slip boundary conditions y = for 9 = 0, a. In 
order to compute (3k, we will argue that the supremum in (13. ip is a maximum and 
use a variational argument. 

The existence of a maximizer is proved by a standard argument in the calculus of 
variations: It is clear that < /3k < oo, and that the ratio I Pt k/I u ,k is a homogeneous 
function of y. Thus we may choose a maximizing sequence of vector functions y with 
fixed H 2 Sobolev norm on [0, a}. Evidently, the quantities L\y\, L 2 y2 remain bounded 
in L 2 , and the complex scalar quantities deyi, dgy 2 at 9 = 0,a remain bounded. We 
may choose a subsequence converging weakly in H 2 such that the quantities deyi, dgy 2 
at 9 = 0, a converge. Then the weak limit is a maximizer by weak lower semicontinuity 
of the L 2 norm. 

Next, consider any smooth curve r i-> y = y(r) into H 2 with the property that 
( 13. 8 p and the no-slip conditions hold for all r, and I p ,k/Iv,k achieves its maximum at 
t = 0. Then at r = we have 

(3.12) = I p> k — (3kl u ,k- 

We now determine I p ^ and I Ut k and solve for /3k- Differentiating I p ^, we find 
' k 2 + a 2 ' 



{a + a + + a + a + + + a_d_)(l — \u\ 



k 

From (13. 8 p we infer that 

a + (l-u J 2 ) = d e (y 1 +y 2 )e^ e '^\l 

(3.13) 



a ^l-u 2 ) = de{ yi +y2)e~^ a) X 



o • 



By differentiating in r, we can solve for a + and d_, allowing us to eliminate a + and 
d_ from the formula for I P) k- Indeed, if we let 7«(6 I ) = doiji for i — 1, 2, we have 

« + = ^§(7i+7 2 )e (fc - )( ^|o, 
1 — u; z r 



^-=lT-^S(7i+T 2 ) e - (fe -^lo- 
|1 — u; z | z 

Similarly, we differentiate J U) fc. Letting L* and denote the formal adjoints of L\ 
and L 2 , respectively, and recalling from the no-slip boundary conditions that y = 
at 9 = 0, cr, we integrate by parts to conclude that 



where 

P = yi( L i L i)Vi + 2/2(^2^2)2/2, ^ = dey^L^i) + d e y 2 (L 2 y 2 )\%. 
We observe that 

L\L\y\ = and L 2 L 2 y 2 = 0, 
so that J^fc reduces to = 4M(if)). 
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Using this information, we can rewrite (I3.12p as 
(3-14) 

= j 7 i (2p k L lVl - (^j^) ^|P(a + e( fc -)( e -) + a^^] 

+ {72 (2/3,L 2 y 2 - [p^-) il^p (a+e ( fc -)^) + a_ e -< fc -*>«) 

Since 71 (6 1 ) and 7 2 (#) are arbitrary at 9 = and 6 = a, (13.141) yields four (natural) 
boundary conditions: 

fk 2 + a 2 \ (l-M 2 ), , J 

2\i k L x y x = I — - — 1 - - - 2 (a+ + a_cu) and 

2p k L 2 y 2 = I — 7 — J 1 _ 2 («+ + a-wj, when = <7, 



£; 2 + a 2 \ (1 



2/3 k Liyx = ( — - — J - _ - 2 (a+u; + a_) and 
2PkL 2 y2 = ( ; — ^(a+ou + a_), when = 0. 



(3.15) 



fe y 1 - w 2 

In addition, we have the four no-slip boundary conditions 
(3.16) yi (a) = y 2 (a) = yi (0) = 2/2(0) = 0. 

3.4. Using fl3~16D . f l3TT5l) . (BSD , and the property L\L\y\ = L* 2 L 2 y 2 = on (0,0-), 
we can explicitly solve for the maximizer of To simplify the calculations in what 
follows, we first use reflection symmetry to show that either 

(a+, a J) = (1, 1) or (a+, a_) = (1, -1). 

Letting 9 = a — 8, we see from our construction of q in (13.41) that a + and a_ exchange 
roles after reflection; thus, it is natural to set a + = a_, and d_ = a + . In addition, 
we let y 2 (6) = yx{6) and y\{9) = y 2 {9). A straightforward calculation shows that 
(yi, y 2 , &+, at-) solves the set of linear equations consisting of ()3.16|) . (I3.15p . (13. 8p . 
and L\L\yx = L* 2 L 2 y 2 = 0. We deduce that 

(yi + Vi, 2/2 + 2/2, «+ + «-, «_ + «+) and (y x y 2 - #2, «+-«-, a_ - a+) 

also solve these equations. We conclude that every pair (a+,a_) will yield the same 
value for (3 k as either (a+,a_) = (1, 1) or (a + ,a_) = (1, —1). Therefore it suffices to 
consider only these cases. 

3.5. We can eliminate y 2 by observing that if (a+, a_) = (1, 1), then y 2 {6) = y 2 {6) = 
y\{a — 9), and if (a + ,a_) = (1,-1), then y 2 {6) = —y 2 {9) = —y\{a — 6). Then we 
infer from boundary conditions in (I3.8P that 

1 - u = d yi(a) - (0), when (a+, a_) = (1, 1), and 
1 + u = 0*2/1 (<r) + 9flj/i(0), when (a + , a_) = (1, -1). 

We are now in a position to solve for 2/1 and ultimately {3 k . We first recall that 
L\ = [ik + 1 — a) 2 + <9g — 2ide — 1 , while, formally, the adjoint of this operator is 
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given by L\ = [ik — 1 + a) 2 + <9 2 — 2idg — 1 . The characteristic polynomials of these 
two operators are 

= (// — (2i — k — ia))(fi — (k + ia)), 
= (/i — (2z + k — ia))(iJL — (—k + ia)). 
Since L\Liyi = on (0,cr), we can conclude that y\{0) takes the form 
(3.18) yi {6) = ai e {k+ia){0 - a) + a 2 e- {k - 2t+ia)e + a 3 e^ k ~ ia)e + a 4 e (k+2i - ta){e - a) 

for some constants a,, 1 < i < 4. The boundary conditions 7/i(cr) = yi(0) = 0, com- 
bined with ( 13. 181) . yield the two equalities 



= ai + c^we 2 " 7 + a 3 co + a 4 , 

(3.19) 

= a\UJ + a 2 + a 3 + a 4 uie m . 

We will use the boundary conditions for y± in (13. 151) combined with the equalities in 
(I3.19P to write the four unknowns Oj, 1 < j < 4, in terms of a + and a_. 
Using the equality Lj = L* + 4(1 — a)ik for j = 1, 2, and (13 .181) . we conclude that 

L Wl = 4(1 - a)ik {a 3 e- {k - ia)e + a^ k+2t - ia) ^ ] ) . 

Plugging this information into the two boundary conditions in (13.151) yields the two 
equalities 

'k 2 + a 2 \ ( 1-M 2 ) 

v .. , 1-u 2 
(3-20) _ 2X , , 2N 



3 fe (l - a)ik(a 3 u + a 4 ) = ( — ) — ^— (a+ + a_a;) 



A; 2 + a 2 \ (l-|w| 



8/3 fc (l - a)ik(a 3 + a 4 we 2mt ) = I - — j I — — ! ^-^(ot_ + a+cA 

\ A; / 1 — ur 

3.6. The value of (3k is determined by the equations in (13.201) and (13.191) together with 
( 13 . 1 TP . Evidently (3k = m&x{(3 + ^, (3-,k} where (3 +t k and /3_^ are the values determined 
from these equations in each of the two cases (a+, a_) = (1,1) and (a+, a_) = (1, —1) 
respectively. 

With (a + ,a_) = (1,1), using the four equations given in (I3.20j) and (13.191) . we 
solve for the unknowns 1 <j < 4, finding that 

a 1 /3 + , fc (l-o; 2 e 2iCT ) = -0 1 (l-a;e 24CT ), 

a 2 /3 + , fc (l-a; 2 e 2iCT ) = -0 1 (l-a;), 

a 3 /3+, fc (l - u;V 2i<J ) = &(1 - a;e- 2 ^), 

a 4 /3 + , fc (l-w 2 e- 2 ^) = 0!(l-a;), 

where 

(k 2 + a 2 )(l-\co\ 2 )(l + co) 
01 8iP(l - a)(l - A) ' 

Using (I3T7D with (I3TT8D . we see that 

1 - w = a x k(\ + a 2 (2i - k)(ue 2i(7 - 1) 

(3.22) 

+ a 3 (ife)(l - Q) + a 4 (2i + k)(l - oje~ 2i(J ), 
where k = k + ia. Plugging the formulas for aj into (I3.22p and solving for (3 + k yields 
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To find /3_ fc, we let (aq_, a_) = (1, —1), and we again use ( I3.20p and (13.191) to solve 
for a,j, 1 < j < 4. To simplify notation, we define 

(fc 2 + a 2 )(l-|q;| 2 )(l-a;) 
8^ 2 (l-a)(l -tu 2 ) ' 
We compute the cij and conclude that 

ai0_, fc (l - w 2 e 2iCT ) = -0 2 (1 + we 2 "), 
a 2 /3_, fc (l-a; 2 e 2 ^)=0 2 (l + a;), 
a 3 /3_ )fc (l - u>V 2iCT ) = -0 2 (1 + W e- 2w ), 
a 4 /3_, fc (l - a;V 2iCT ) = 2 (l + a>). 
We solve for using (13 . 1 7j) with (13 . 18[) like before, and find that 

^' (l-w 2 e 2iCT ) V ; (l-d; 2 e- 2lCT ) V ; J 

At this point, one can check that (3±^ = j}±^ as given in (11.9j) . 



3.7. To complete the proof of Theorem ll.4[ as indicated at the beginning of this 
section, we must show that f3 aj0t > f} a>a . To prove this, suppose (3 < $ a>a . Then there 
exists k ^ such that /3 fco > (3. We choose y to be a maximizer of the ratio I Pl k /I u ,k - 
In a change of notation, we let I P) k and J Ui fc denote the integrals corresponding to this 
fixed y, with q determined by (j2.9j) for varying. Since y may not be a maximizer for 

7^ feo, we only have /3& > I P} k/Iu,k m general. However, by continuity it is evident 
that there exists 5 > such that whenever \k — k \ < 5 we have I p ^/I u ,k > $■ 

Next, we define x<s(^) to be a smooth bump function independent of 9 and supported 
in a 5- neighborhood of fco- Recalling that — iw = Vy, we set ws = Xsib an d Qs = X<5<?> 
and we observe that (wg, qs) solves f!2.9p and ws = for 9 = 0, a. Moreover, if I Ps ^ and 
I Ugt k are the integrals corresponding to ws and qs, then one sees that I PSt k = Xs 2 I P ,k 
and I USt k = Xs 2 Iu,k- We can then write 

(3.24) '» tt-jy'** dk - ; Csj^* dk j. 

J ^ //J^ X<5 2 4,fc ^ /^J"/ X6 2 Iu,k dk 

We conclude that /3 CT>a > /3, hence /3 CT)Q > /3 a%a . 

4. Causes for blowup of the optimal constant 

One can rewrite the formulas for /3 + & and f3_ ^ from Theorem 11.41 in the following 

way: 

(A l) a = ^1+^2 

1 ' } P± ' h 2k 2 (cosh(ka) =f cos(aa))(cosh(2A;a) - cos(2(l - a)a)) ' 

where 

= (k 2 + a 2 ) sinh(fco") [sinh(2A;<7) =F 2 smh(ka) cos(a) cos((l — ot)a)\ 

and 

k(k 2 + a 2 ) sinh(/cer) r / , N s , s tt , 

V> 2 = — ^ — - sin(2(l-a)o-)=F2cosh(fccr) sin((l - aW) cos(cr) . 

1 — a 1 
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From (14. ip it is clear that for fixed a ^ 1 and fixed a G (0, 27r), and p-^ as 
functions of k are continuous everywhere except k = 0. If we take the limit of (14.11) 
as k approaches 0, we find that 

(4.2) li.n i:., ?±,0 



where 
(4.3) 



k^o ' ' 2(1 =f cos(aa))(l - cos2(l - a) a) 



ip3 = 2a =F o (cos(aa) + cos(2 — a)a) 



a 



1 — a 



^sin(2 — 1a)o =p (sin(2 — a)o — sin(acx))). 



From (14.21) we see that P aa typically blows up when either aa = mr or (1 — a)a = mr 
for some n G Z. 

The first set of singularities above is a result of the unboundedness of the Neumann 
problem for the Laplace operator in weighted spaces on a cone. To see this, we observe 
that in (I3.13p . a + and a_ become undefined as k — > when u 2 = e ~ 2( ~ k + ia ) u — y i ( 
which occurs precisely when aa = mr for nGZ. 

The second set of singularities above, which occur when (1 — a)a = mr, result 
from failure to bound the boundary data n ■ (A — VV-)w in terms of Am. For these 
combinations of a and a, the L 2 norm of r a Au in JC a is not sufficient to control 
n ■ (A — VV-)m appropriately. This is fundamentally due to the existence of harmonic 
fields u = cr x ~ a sin((l — a)9) where c is a constant vector. Corresponding to these 
fields, there are nontrivial modes (2/1,2/2) for k = satisfying L\yi = = L22/2 and 
the no-slip boundary conditions (I3.16p . while (a + ,a_) is non-zero. One finds then 
that the maximum of I P) k/Iu,k ~ > 00 as k — > 0. 

To see just how this can occur in terms of the computations of section 3 for certain 
combinations of a and a (7^ 1 or 0), we observe that in section 3.5, Lij/i = iff yi{0) 
takes the form given in (I3.18P with = a 4 = 0. One can then satisfy the no-slip 
boundary conditions through (13.191) for some nonzero a±, a 2 if and only if u 2 e 2la = 1, 
meaning k = and (1 — a)a = mr for some n G Z. We may simply take 2/2 — 0, and 
it follows by (13. 3 p that I Ut o = 0. 

But then, a\ = — g^w 7^ 0, and we compute that dgyi = 2i(a — l)a± 7^ at 9 = a, 
yielding nonzero boundary values for dgq in (13. 7p and causing I p ^ to be positive in 
(12. 3p . If we vary k while holding (2/1,2/2) fixed and use (13. 8 p to determine (a + ,a_) 
and thence q, we see that — » as k — > while —> / p ,o > 0. This results in 
j3 k — > 00 as k — > 0, hence {3 a>a = 00 when (1 — a)a = mr. 



5. Proof of Theorem 11.21 

We now use Theorem 11.41 to prove Theorem 11.21 through a localization argument. 
Let fl denote a bounded domain with a straight corner. Replacing Q by a suitable 
rotated translate if necessary, we may assume there is a neighborhood U of such 
that U fl O = U H /C CT , where 7^ ir. 

Fix any < 1 and CeR. We observe from the formula for /3± £ with a = given 
in (11.1 Op that /^o > 1 when a 7^ ir . Therefore, there exists a solution (u,p) to (12.81) 
with u in C^ c (/C (7 \{0}, M 2 ) which satisfies J K |Vp| 2 > /3 J K \Au\ 2 . Replacing (u,p) 
by suitable dilates if necessary, we may assume that the support of u is contained in 
U. 
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We construct a sequence of solutions (uj,Pj) to (12 .8p on Q by setting 

Uj (x) = r X u{jx)\u and Vpj = (I - P)(A - VV-)u,- in O. 

We see that Apj = in Q and n ■ Vpj = n ■ (A — W-)uj on dQ. Moreover, since 
Uj is supported in fl fl K a for all j, we have || Ai^H^^) = || Aw|| L 2(jq) < || Aw]!^^) 
and || Vttj||z,2( n ) = j -1 1| VttH^^) < .7 _1 || VttH^^) for every j. This construction 
allows us to write the following series of inequalities for sufficiently small e > and 
for sufficiently large j: 



[ {(3\A Uj \ 2 + C\V Uj \ 2 ) < [ ((3\Au\ 2 + Cj- 2 \Vu 

JO, JKa 



2^ 



< / (P\Au\ 2 ) + e< / \Vp\ 2 

We claim that 

|2 



(5.1) / |Vp| 2 < liminf / |Vp. 

JK a Jn 

To see that ( 15.1ft holds, we first use the definition of Uj, the equality Vpj = (I 
P)(A — W-)uj, and orthogonality of the Leray projection to observe that 

(5.2) ! \V Pj \ 2 < [ \(A-VV-)u\ 2 < [ |(A-VV-)m| 2 
Jn J jn JK a 

for all j. If we define 

I x\ I f fx 
( 5 - 3 ) P*A X ) =Pi 7 ~ 7775T / Pi 



jj m(B) J B \j 

for x G j^l, where B corresponds to the domain B given in (jl.7|) . then we can apply 
a generalized Poincare Inequality (see [51 Ch. 2]) to conclude that for each n G N, 



(5-4) / \ P *\ 2 <C n / \Vp*\ 2 <C n / \V Pj \ 2 

Jnn Jnn Jn 

for sufficiently large j. By a standard diagonalization argument, we can construct 

a subsequence of {Vp*}, which we henceforth denote as {Vp*}, converging weakly 

in L 2 (nfl) for every n G N. This implies by (15.4ft and by another diagonalization 

argument that, up to subsequences, {p*} converges weakly to some p* in L 2 (nfl) for 

all n G N. By uniqueness of weak limits, we can conclude that {Vp*} converges 

weakly to Vp* in L 2 (nQ) for every n. Moreover, by properties of weakly convergent 

sequences we can write f nU |Vp*| 2 < liminf j^,^ f nn |Vp*| 2 for each n. We can then 

conclude that for sufficiently large j, 



Kn 



|Vp*| 2 < lim / |Vp*| 2 < liminf / |Vp*| 2 = liminf / \Vpj\ 2 . 
n ^°°Jnn Jin 3-+°° Jn 



It remains to show that J K |Vp*| 2 = J K |Vp| 2 . This will imply ( 15. lft . 

To show that j K |Vp*| 2 = L. | Vp| 2 , we first show Ap* = in K, a . We fix a compact 
subset K of /C CT , and we apply the mean value property and weak convergence of {p*} 
to conclude that for any y G K, \pj(y)\ < C\\p* ||L 2 (n^) < C, giving equiboundedness 
of {p*} on K. Moreover, by the mean value property and weak convergence of {Vp*}, 
{Vp*} is equibounded on K, implying that {p*} is also equicontinuous. Therefore, 
up to subsequences, {p*} converges uniformly on K to p*. We again apply the mean 
value property and uniform convergence of {p*} on K to conclude that p* is harmonic. 
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Since Ap* = on nQ, we infer that the sequence {Vp*} converges weakly to Vp* 
in H(div,nQ), the space of vector fields in L 2 (nQ) with divergence in L 2 (nQ). By the 
boundedness of the trace operator mapping H{div, nQ) into H~2(d(nQ)) (see, for 
example, [I] Theorem 2.5), we can conclude that n ■ Vp* converges weakly to n ■ Vp* 

in H~^(d(nQ)). As n ■ Vp* = n ■ Vp on dfC a n d(nQ) for every n, it follows that 
n ■ Vp* = n ■ Vp on d)C a . 

Using the equalities Ap = Ap* = in K a and n ■ V(p* — p) = on dJC a , we can 
now integrate by parts to conclude that J K |V(p* — p)\ 2 = 0. For e C^°(/C cr ), we 
have that 

/ V<f)-V{p* -p) = <pn-V{p* -p) - <pA(p* - p) = 0. 

Since p*—p belongs to Y and C^°(JC a ) is dense in Y, it follows that j K \ V(p*— p) | 2 = 0, 
and ( 15.1 1) holds. This completes the proof of Theorem 11.21 
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